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Abstract 

We study the constraints imposed by conformal symmetry on the equations of fluid dynamics 
at second order in gradients of the hydrodynamic variables. At zeroth order conformal symmetry 
implies a constraint on the equation of state, £q = |P, where £q is the energy density and P is the 
pressure. At first order, conformal symmetry implies that the bulk viscosity must vanish. We show 
that at second order conformal invariance requires that two-derivative terms in the stress tensor 
must be traceless, and that it determines the relaxation of dissipative stresses to the Navier-Stokes 
form. We verify these results by solving the Boltzmann equation at second order in the gradient 
expansion. We find that only a subset of the terms allowed by conformal symmetry appear. 



I. INTRODUCTION 



There has been great interest recently in the dynamics of conformally invariant fluids. 
One motivation is the observation that the AdS/CFT correspondence [1], which relates the 
dynamics of a four-dimensional conformal field theory to quantum gravity in ten dimensions, 
also implies a fluid-gravity correspondence [2-4]. The fluid-gravity correspondence relates 
solutions of the Navier-Stokes equation for conformally invariant fluids in four dimensions 
to dynamical black hole solutions of the Einstein equations in higher dimensions. The fluid- 
gravity correspondence has been used to establish the form of second order terms in the 
equations of relativistic fluid dynamics [3, 5]. The correspondence has also been extended to 
certain non-conformal fluids [6], superfluids [7, 8] , and to the incompressible non-relativistic 
Navier-Stokes equation [9]. 

Another motivation comes from the experimental discovery of nearly perfect fluidity in 
certain conformal or almost conformal fluids [10]. The first example is the observation of 
nearly ideal flow of the quark gluon plasma produced in heavy ion collisions at the relativistic 
heavy ion collider (RHIC) [11]. This discovery was soon followed by the discovery of nearly 
ideal fluid dynamics in a very different system, the dilute Fermi gas at unitarity [12]. Both 
the quark gluon plasma and the dilute Fermi gas are characterized by a very small shear 
viscosity to entropy density ratio, i]/s < OAh/ks [13-19]. This result is close to the value 
rj/s — H/(4:7ikB) obtained in the strong coupling limit of a super-conformal Yang Mills 
plasma [20, 21]. 

In the present work we will focus on the dilute Fermi gas and study the hydrodynamics 
of a scale invariant non-relativistic fluid. At unitarity the two-body scattering length of the 
fermions is infinite and the range of the interaction is zero. This implies that the system 
is exactly scale invariant [22-24]. Hydrodynamics is an effective theory that governs the 
behavior of the system at long distances. The equations of motion are organized as an 
expansion in the number of derivatives acting on the hydrodynamic variables. In the case of 
a one-component non-relativistic fluid the hydrodynamic variables are the temperature T, 

— * 

the density n, and the fluid velocity V. At zeroth order in the gradient expansion conformal 
symmetry implies a constraint on the equation of state, So = |P, where £o(n,T) is the 
energy density and P(n, T) is the pressure. At first order, conformal symmetry implies that 
the contribution of the divergence of the velocity field to the pressure, 5P = £(V • V), must 



1 



vanish [25-27] . This means that the bulk viscosity ( of a scale invariant fluid is zero. In this 
work we will study constraints from conformality at second order in the gradient expansion. 

Second order hydrodynamics has a long history [28-30] but the theory is rarely used in 
practice. An important exception is the relativistic case, where the second order theory has 
the advantage of being manifestly causal [31]. There is a similar reason for using second order 
hydrodynamics in certain non-relativistic flows. The expansion of a scale invariant fluid after 
the release from a harmonic trap leads to a Hubble-like flow in which the expansion velocity 
is linear in the coordinates [18]. This implies that ideal stresses propagate outward with a 
finite velocity, but at first order in gradients dissipative stresses build up instantaneously 
everywhere in space. This unphysical behavior can be avoided by taking into account a 
finite relaxation time for the dissipative stresses. We will show that relaxation time effects 
appear naturally at second order in the gradient expansion. 

There are additional reasons for studying higher order terms in the equations of fluid 
dynamics. Higher order terms provide information about the convergence of the gradient 
expansion, and they lead to new Kubo relations. In the relativistic case it is known that 
some of these Kubo formulas relate transport properties to thermodynamic quantities [32]. 
Finally, understanding the structure of the gradient expansion may help in constructing 
holographic duals of non-relativistic conformal field theories [33-35]. The paper is structured 
as follows. In Sect. II we review a local extension of Galilean and conformal invariance and 
study constraints on the structure of gradient terms in fluid dynamics. In Sect. Ill we 
consider a specific kinetic equation and solve for the stress tensor and the entropy current 
at second order in gradients. We show that the result agrees with the general constraints, 
but that the kinetic theory only generates a subset of the allowed terms. We conclude with 
a discussion of open questions. 

II. CONFORMAL HYDRODYNAMICS 
A. Navier-Stokes equation 

The constraints imposed by Galilean invariance and conformal symmetry on the structure 
of the Navier-Stokes equation were studied by Son [25] based on earlier work by Son 
and Wingate [22]. The basic strategy is to generalize Galilean and scale invariance to 
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local symmetries. For this purpose we consider the fluid moving in a curved background 
characterized by the metric gij(x,t). We also include background U(l) gauge fields A (x,t) 
and Ai(x,t). The Navier-Stokes equations in curved space are 

^=dt (y/gp) + V, (pV*) =0, (1) 

^pd t (y/gpVJ + V k U\ = —Qi, (2) 

^jft (y/gs) + V kJ k s = f , (3) 

where g = det(^j) and Vj is the covariant derivative associated with the metric gij. The 
hydrodynamic variables are the mass density p, the entropy density s, and the fluid velocity 
V 1 . The energy momentum tensor is 11^ = IT°- + SILij, where 

n° = pViV 3 + Pg l3 (4) 

is the ideal fluid part, and <SIIy is the viscous correction. The entropy current is f s = sV % +5f 8 , 
and R is the dissipative function. We will specify the viscous terms 6Hij,5f s and R 
below. The Lorentz force Qi is defined as Qi = Ei — F^V^ with Ei = d t Ai — diA and 
E{j = diAj djAj,. 



B. Diffeomorphism invariance 



The equations of fluid dynamics are invariant under a combination of local 
diffeomorphisms x % — > x % + t) and gauge transformations A — > A — a, Ai — >■ — c^a. 
This symmetry generalizes Galilean invariance. A Galilei transformation with boost velocity 
U l corresponds to a flat background metric g^ = Sij and £ l = C/ l i as well as a = mU % x l . The 
hydrodynamic variables are invariant under gauge transformations. They transform under 
diffeomorphisms as 

5 P = -e fe V feP , (5) 
5s = -d k V k s, (6) 
SV l = -^VhV* + V k V k C + C . (7) 

We note that p and s transform as scalars, but Vj does not transform as a vector. We have 
5V l = C^V 1 + £\ where £ ? is the Lie derivative and £ l is an "extra" term. The gauge fields 
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and the metric transform as 



5A = -eV k A - A k i k , (8) 

SAi = -C fe V fe A - A k V t e + mgd k , (9) 

Sga = -QikVj^ -g kj V£ k . (10) 

Using the transformation laws for the gauge fields we can deduce the transformation 
properties of the field strengths and the Lorentz force. We find 

5E t = C 6 E Z + F lk i k + mg lk i k + mg tk £ k , (11) 

5F tj = C^j + mg jk Vii k - mg ik V £ k , (12) 

5Q, = C^Q t + m (gd k + 9 lk ^ + V k Vd h - 9i k V j V . (13) 



Using these ingredients we can verify that the hydrodynamic equations (1-3) are invariant. 
We note that under diffeomorphisms the ideal parts of the momentum and entropy currents 
do not transform as vectors, and the ideal part of the stress tensor does not transform as a 
tensor. The extra terms in the transformation laws are canceled by terms that arise from 
the time derivatives of the conserved charges and the Lorentz force. In order to maintain 
diffeomorphism invariance beyond the level of ideal hydrodynamics dissipative corrections 
to the entropy current and the stress tensor have to transform as vectors and second rank 
tensors, respectively. 

At first order in gradients dissipative terms involve the heat current and the shear and 
bulk stresses [36]. The heat current 

q i = -V i \og(T), (14) 

where Vj is the covariant derivative, transforms as a vector under diffeomorphisms. The 
shear and bulk stresses can be promoted to second rank tensors by adding extra terms 
involving gij [25]. We will show below that these terms automatically appear in kinetic 
theory. We have 

_ 2 

<Tij = ViVj + VjVi + gij - -gij{cr) , (15) 
(a) = V-V+f g (16) 

where is the shear stress tensor and gij(cr) is the bulk stress tensor. 
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C. Conformal invariance 



Conformal transformations are infinitesimal rescalings of the time variable, t — > t+f3(t). A 
non-relativistic scale transformation corresponds to j3 = bt combined with a diffeomorphism 
generated by £ fc = ^bx k . Special conformal transformations are generated by (5 = —ct 2 , 
£ fc = —ctx h and a = —\mcx 2 . An operator O is said to have conformal dimension Ao = [O] 
if 

SO = -(30 - Ao/30 . (17) 
Examples of operators with well-defined conformal dimension are 

[g tJ ] = -l, [g ij ] = +l, [A ] = +l, [A]=0. (18) 

We also have [0]0 2 ] = [Oi] + [0 2 ] and [diO] = [O], but d t O does not, in general, have a 
well-defined conformal dimension. The conformal dimensions of the hydrodynamic variables 
are 

w = |, w = |, m=o. as) 

We can now check under what conditions the equations of hydrodynamics are conformally 
invariant. Two key relations are 

Sgij = -Ma + kij = S A=0 (g tj ) + ha , (20) 

and 

where 5 A ~ d O denotes the transformation of a conformal operator O of dimension d. Using 
equ. (20) and (21) we can show that the equations of fluid dynamics are conformally invariant 
if 

M = |. [(J.)i] = l> [R] = l- (22) 
These relations are satisfied in the case of ideal fluid dynamics (in ideal hydrodynamics 
R = 0). At one derivative order we find [a^] = 0, [<&] = and 5(a) = 5 A=1 (a) + \ This 
implies that 

5Uij = -rpij , 5(j s )i = Kqi , (23) 

are conformally invariant provided [rf\ = [k] = |, but 5Uij = —(gij(a) violates conformal 
invariance unless ( — 0. We note that the conformal dimension of r\ and k agree with their 
naive scaling dimension. 
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D. Second order terms: Building blocks 



In order to construct the most general second order terms we begin by listing conformally 
invariant building blocks, that is diffeomorphism invariant scalars, vectors, and tensors with 
well-defined conformal dimensions. These objects are constructed from the hydrodynamic 
variables T, P and V\ the metric and the gauge fields A , Ai, and time or spatial 
derivatives. Instead of T and P it is sometimes useful to consider T and the dimensionless 
variable a = /i/T, where /x is the chemical potential, as independent quantities. We begin 
with one-derivative scalars. Diffeomorphism invariant one-derivative scalars constructed 
from the hydrodynamic variables are DT, DP and (a), where D = d + V k Vk is the 
comoving derivative. Because there are two scalar equations of motion only one of these 
objects is linearly independent. Neither one of the one-derivative scalars has a well-defined 
conformal dimension but we can construct conformal scalars by taking linear combinations. 
The quantity 

S=(D + l(a) > )T (24) 

is a conformal scalar with dimension [S] =2. We will show, however, that at leading order 
in the derivative expansion S vanishes by the equations of motion. We can construct three 
conformal vectors 

V} = V t T, V? = V t P, V? = DV l -^, (25) 

m 

with weights [V/] = 1, [Vf] = §, and [Vf] = 1. There is one vector equation of motion and, 
up to higher order terms in the derivative expansion, V 3 can be expressed in terms of V 2 . 
There are two one-derivative tensors. The first is the shear tensor defined in equ. (15). The 
second is related to the vorticity of the fluid. We have 

73 = <7y, 7j = % = V^-V^-^, (26) 

with [73] = [7- 2 ] = 0. 

E. Second order terms: Stress tensor 

In this section we list conformal two- derivative tensors that can contribute to the stress 
tensor in a conformal theory. Rotational invariance implies that the stress tensor has to be 
symmetric. We will argue below that in a conformal theory the dissipative part of the stress 
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tensor also has to be traceless. One class of two- derivative tensors arises from contracting 
two one-derivative tensors. We find 

0]j = a (i k(T j)k , 0% = a (i k ^j)k , 0% = tt(i k ^j)k • ( 27 ) 

with [Ojj] = \0%\ = [Ofj] = 1. Here, (ij) = \{O i3 - + 3i - \g i3 O k k ) denotes the symmetric 
traceless part of O i3 . There is a unique diffeomorphism invariant tensor of conformal 
dimension +1 that can be constructed from the comoving derivative of a i3 . This tensor 
is given by 

0% = g ik a) + V k V k a l3 + | {a)a i3 + F {l k a kj) , (28) 

where O^j) = \{O i3 + Oji). There is a set of conformal tensors that can be obtained as 
the tensor product of two one-derivative tensors, or as two covariant derivatives acting on a 
scalar. We find 

0% = V (l TV,)T , 0% = V (l PV,)P , Oj 3 = V { {TV 3) P , (29) 
Of.^V^T, 0% = V { N 3) P, (30) 

with [Cfj -9 ] = {2,5, |,1,§}. Finally, in curved space there is a two-derivative tensor 
constructed only from the metric. This is the traceless part of the Ricci tensor, 

0$ = R {ij) . (31) 

We can compare this list with the analogous result in the relativistic case [5]. In relativistic 
fluid dynamics there are five second order terms (four in flat space). Three of those are 
analogous to Oj, — Of,, and one is similar to Of,. It is interesting to note that the relativistic 
analog of equ. (28) has a coefficient | in front of (a) instead of |. This difference is similar 
to the difference between the relativistic and non-relativistic conformal equation of state, 
£ = |P versus £ = \P. In a relativistic theory without conserved charges there are 
only two hydrodynamic variables, T and the four- velocity C/ M , and as a result there are no 
terms analogous to Ofj, Oj 3 and Ofj. Finally, in the relativistic case there is only one linear 
combination of Of 3 and Of 3 which is conformal, and that linear combination can be related 
to Ofj using the equations of motion. We can also compare our result to the structure of the 
Burnett equations [30]. The Burnett equations contain six independent two-derivative terms 
in Iljj. The reason that the number of independent terms is smaller despite the fact that 
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there are no restrictions from conformal symmetry is that the Burnett equations are derived 
from kinetic theory, which only leads to a subset of the terms allowed by the symmetries. 

We end this section by showing that there are no trace terms in SILij, despite the fact 
that second order terms like gij<Jki<J kl are conformal tensors. In Sec. Ill A we will show that 
in a non-trivial background the equation of energy conservation is given be 

±=dt {^gS) + V k3 k e = -\ 9i ^ , (32) 

where 8 = S + \pVjV % is the energy density and j\ = (£ + P)V k + 5j k is the energy current. 
Equ. (32) is conformally invariant if 

£ = |p, g l 'SU t3 =0, (33) 

which shows that the equation of state has to be conformal, and the dissipative part of the 
stress tensor must be traceless. 



F. Second order terms: Entropy current 

In this section we will use the results of Sec. II D to construct two- derivative conformal 
vectors that can contribute to the entropy current. There is one vectors that can be written 
as the comoving derivative of a one-derivative vector 

Ql=Dq t , (34) 

with [Qj] = 1. Note that DViCt is not an independent vector, because it can be related to 
Q\ and the vectors Qf 5 defined below using the hydrodynamic equation for Da. We find 
four vectors that can be written as contractions of one-derivative vectors with one- derivative 
tensors 

Qi = , Q- = Q^ij , Qt = (VaK , Q\ = (Va)^- , (35) 

with [Q^~ 5 ] = 1. Finally, there are two vectors that can be constructed from the covariant 
derivative V« acting on one-derivative tensors, 

shv^, Qj = vm tJ , (36) 

with [Q^ 7 ] = 1. Again, the number of terms is bigger than the number of independent 
parameters in the entropy current in the Burnett equation, which is five. 
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III. KINETIC THEORY 



A. Conservation laws 

In this section we will investigate the structure of higher order terms starting from the 
Boltzmann equation. Our motivation is three-fold: i) We will verify that the results indeed 
satisfy the constraints derived in the previous section; ii) we will show that kinetic theory 
(within the relaxation time approximation) only leads to a subset of the allowed terms; iii) 
we will demonstrate that this subset manifestly satisfies the second law of thermodynamics. 

The kinetic equation in a curved background g^ (x, t) can be found by starting from the 
Boltzmann equation in a four-dimensional curved space [37-39] 

jo (V^r - r >v^) f<t>*>p> = °w > ( 37 ) 

where f(t,x,p) is the distribution function, C[f] is the collision term, F" is the Christoffel 
symbol associated with the four-dimensional covariant derivative V M , i,j, k are three- 
dimensional indices and /i,a, (3 are four-dimensional indices. In the non-relativistic limit 
[22] 

r>V = / - -£ W + -g il F lk p k p° + g a g lkP k p° , (38) 

where p° ~ m, E l is the electric field and Fi k is the magnetic field introduced in Sect. II A. 
In the following we will drop the Lorentz force (these terms can always be restored using 
the symmetries described in Sect. II B). We get 

Vf(t,x,p) = C[f], (39) 

where we have defined the Boltzmann operator 

at mox k m dp 1 op 1 

Consider a collision term that describes elastic two-body collisions p\ + p 2 — > p?, + p±. The 
symmetries of the collision term imply that 

J dF^ X C[f] = l -j dF^fC[ X i + X2- Xv - X2>] , (41) 

where dr = d 3 p and x% = xipd- This implies that moments of the collision term with regard 
to the collision invariants x — m i P\ 9%jV % V^ I (2m) must vanish. We will denote averages over 



the momentum distribution as 

(A} = JdT^gfA. (42) 
We define the conserved charges and currents as 

P=(m), n k = {p k ), g = -L(g tjP y), (43) 

= ^ W ) , Jl = <pW/> ■ (44) 

The three collision invariants lead to three conservation laws. The first conservation law is 
obtained by taking a moment of the Boltzmann with x — m. We get 



= J dr y/gmVf(t,x,p) 



{m) l(m) + d k </> - |f </> - A <i>y + m,W> + (21V + ^ 



^ 2 (? v ' ' / 2(7 x ^ 7 dp* x ' » »™ 1 ' \ ' g 

The momentum equation follows from taking a moment with x — dijP 3 



= | dr x>/(* , x, P ) = ^ + v*n* . 



(46) 



° = 2^ / ^ V^^V' 2>/(t, x, p) = ^ j + V fc j* + ^TF . (47) 



and the energy follows from integrating over x — dijP 1 ^ /(2m), 

This is the result we have used in equ. (32) above. In order to determine the form of the 
currents in terms of the hydrodynamic variables we need to know the functional form of the 
distribution function near equilibrium. We will see shortly that f(x,p, t) = f^^g^dc?) with 
Ci(x, t) — Vi — Vi(x, t) and Vi — Pi/m. Then 

vr, = pV, , IUj = pViVj + P 9ij , i =(e + P + \pVjV^j V 1 , (48) 

with S = \(m(v — V) 2 ) and P = |£ - We can now use the thermodynamic relation 
d£ = Tds + ^dp to determine the entropy equation. We get 

M^f) +Vfc ( a y*) = o. (49) 

This result can also be derived directly from the Boltzmann equation. Computing 
the moment of the Boltzmann equation with regard to log(/) (for classical particles, 

10 



log(//(l ± /)) for Bose/Fermi statistics) we obtain equ. (3) with the entropy density given 
by 

s = J dT^/g 
and the dissipative function 

f = / dr ^i°«(T^7)c[/], (51) 

with a = 0, ±1 for classical, Bose and Fermi statistics, respectively. In the case of 
classical statistics equ. (51) implies that R vanishes provided the distribution function is 
an exponential of the collision invariants. The same result holds for quantum statistics, 
with the exponential replaced by Bose-Einstein or Fermi-Dirac distribution functions. This 
implies that 

f<°Ht,x,p) = /<°> ( ^^^ t,X) - «(*,*)) , (52) 
where is the Maxwell, Bose-Einstein, or Fermi-Dirac distribution and a = pi/T. 



/log 



l + af 
f 



+ a log (1 + af) 



(50) 



B. First order solution 

In this section and the next we will obtain a solution of the Boltzmann equation to second 
order in the gradients of the hydrodynamic variables. We write the distribution function as 

/ = / {0) + Sf = /<°> + /W + + . . . , (53) 

where /W contains terms with i gradients. Since we are interested in the structure of the 
terms that appear in 5f, and not in computing the values of transport coefficients for a 
specific theory, we will consider a very simple choice for the collision term, the relaxation 
time approximation 

C[f] = - 6 -l- (54) 

Conformal invariance implies that A must have conformal dimension 1. This is physically 
reasonable, since A has units of time. The collision term conserves mass, momentum, and 
energy. These conservation laws constrain the form of Sf. We have 

Jd Xc Sf = 0, Jd X Jfd = 0, J d Xc Sf Qijc'cP = , (55) 
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where we have defined dx c = \/# d 3 c/(2-7r) 3 . The constraints imply that the correction to 
the stress tensor due to 5f is given by 



511*' = J & Xc 5fpy = J ■ 



d\ t .6fi>'i>> = I d\,.Sf (cV'-^V 



(56) 



where c 2 = <?yCV and from now on we will use units m = 1. In order to compute 
we will follow the procedure of Chapman and Enskog [40, 41] and compute the LHS of the 
Boltzmann equation using the zero'th order solution equ. (52). We find 



vf(0) = - /( ° )(1 2 + a/(0)) ^ + 2&cPg ij +p k c i cPd k g ij + 2p k ePg ij & l 



Using di = 
vf (o) 



- 2g lk T) lT Pp l c k - 2g kj p*c k - <^cV log T - p fc g fc ) - 2T (d + p fc 3 fc a) ) . (57) 

-d t Vi, djCi = —djVi, and p i — a + Vj we obtain 
/<°>(l + a/<°>) 



2T 



(2Td + 2TU fc <9 fc a) + 2d (V* + U fc V fc U, + Td.a - ^) 



+ c z ( $ logT - U*g, + - ( V • V + 



+ 



3 y 



2^ 

2 / q 
V 3 V k + V fc V) + g kj - -g jk f V • V + ^ 



c 2 c fc g fc + aic fc g fc > , 



(58) 



where, in order to make the constraints equ. (55) manifest, we have added and subtracted 
the term a\c k q k . The first three terms in equ. (58) are proportional to the collision invariants 
and must vanish by the equations of motion in order to satisfy the constraints. The fourth 
term, which is proportional to QCj — ^c 2 gi k , automatically satisfies the constraints. The fifth 
term, proportional to c 2 c k , satisfies the constraints when combined with a 1 c k term. This 
implies 

|d Xc / (0) (l + «/ (0) ) (c 4 -a lC 2 )=0. (59) 

Using the thermodynamic identities derived in App. A we find ai(t,x) = 5P(t,x)/n(t,x). 
For a Maxwell gas the equation of state is P = nT and a\ = 5T. The terms proportional 
to the collision invariants 1, q and c 2 involve time derivatives of the hydrodynamic variables 
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a, T and V*. They are easily seen to vanish by the Euler equations 

= n + V k d k n + n (v ■ V + J^J , (60) 

= Vi + V k V k V t + ^ , (61) 

n 

= g + V h d k S + (£ + P) (V • V + ^ . (62) 

We can now solve the first order Boltzmann equation X>/(°) = C[/(°) + /^] for Using 
equ. (54) we find 

(1) A/(°)(l + a/W) , < , , 2 , fc , 
r j = 2T ^ * j ~ ai ' ° qk > ' ^ ' 

where is the shear tensor defined in equ. (15). Using the result for f^> we can compute 

the dissipative corrections to liij and f s . We find 511^ = —r)(Jij and 5f s = nq l . We note that 

the solution of the Boltzmann equation in curved space automatically leads to a conformal 

and diffeomorphism invariant stress tensor and entropy current. The shear viscosity and 

thermal conductivity are given by 

A 



rj 



K = 



15T 
A 



|d Xc / (0) (l + «/ (0) )c 4 , (64) 
12T2 -/dx c / (0) (l + a/ (0) )(c 2 -«i)c 4 . (65) 

Using the results in App. A we can show that rj = \P and k = -^(7Q — 75P 2 /n), where 
Q = (c 4 ). In the case of a Maxwell gas the result for the shear viscosity reduces to the 
familiar form 77 = XnT. In this limit we also find k = |AnT, which corresponds to a Prandtl 
ratio Pr = c p r}/n = 1. Note that this result is a consequence of the simple form of the 
collision term in equ. (54). In a more complete treatment the high temperature limit of the 
Prandtl ratio is Pr = | [42]. Finally, the dissipative function R is given by 

R = + K Tg i q i , (66) 

which is manifestly positive for A > 0. 

C. Second order solution 



At second order in the derivative expansion the Boltzmann equation reduces to 

Ki) 1 ( r>f(°) • ^ (2) 



Vf^+[vr + J —\=- J —. (67) 
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We note that (Z>/(°) + jf^) vanishes at first order. At second order we have to include 
gradient terms in the equation of motion and this term is not zero. The Boltzmann operator 
acting on gives 

A/<°>(l + a/<°>) f„, , , , , ,2 



2T 



-j/3(cV<7y- (c 2 -ai) c k q k )' 



+ (cVay - (c 2 - ai) c fe g fc ) (z> log(A) - D (logT) ) 
+ V (c i c'(r lJ - (c 2 - ai) c fc g fc ) | , 



(68) 



where we have defined /3 = (1 + 2af^)/(2T). We will decompose this expression as 

P/« = (P/) ortft + (Vf) scal + (Vf) zm , (69) 

where (T>f) ort h are terms that are orthogonal to the zero modes of the collision operator, 
{Vf) sca i is orthogonal to the zero modes but does not contribute to the conserved currents 
IT- 7 and f e , and (Vf) zm contains terms that are proportional to the zero modes. This means 
that (Vf) zm must cancel against (X>/(°) + jf^)- Before we check this we collect the terms 
in (Vf) orth . We get 



(ww=- A/(0,(1 „; a/ '°V 



2T 



C t C j C k C l O r ij(Ji c l — — C 4 (T 



15 



c l c j q i q j 



c 2 q 2 



+ 



c A c l c 3 q i q j - 



c 6 q 2 



- 2ai/3 



c 2 c t c j q i q j - 



c 4 q 2 



(3c c l c 3 c Oijq k —cq k a, 

o 



+ 2ai [/3c i c j c k a ij q k - c l q k a k ] - 2 
+ [cV<7 - - c 2 c l qi + a^qi] (D log(A) - d t logT + V k q k ) 



2 ' . L . J.. ~k 



+ 
+ 



i 4 k 2d\ ■ j, 

c l c J c k a ij q k - —c l q k a k 



cVc fc _ 2oi c , jk 
5 



c 4 g 2 



+ 



a 1 c l c j q i q j 



a\C 2 q 2 



+ [cV - aic fe ] ^VfcFi + W - q k - V l V t q k + 



c 2 cV - 

2(a) 



Viqj 



Qk 



+ 



cV - 



9 ij c 2 



g jk a k + V k V k a i:j + a^qj - a ik a) - Q ik a* - 



2a ij {a) 



+ 2 [Vi + V K V k Vi qj + fcdj-ai 



(70) 



The terms in the square brackets are automatically orthogonal to the zero modes. As in 
the previous section this is achieved by adding and subtracting terms. Subtraction terms 
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that are proportional to the zero modes are collected in (Vf) zm . The remaining terms are 

collected in (Vf) sca i. We have 

A/<°>(l + a/<°>) (20o* 2 3a\q* 2 2 , 

■S -jg- [c - 6 4 ic - O40J H — [c - o 2 iC - 620J 

2(3a iq 2 



scaZ 



2T 



&61C 



&60] 



641 c 



3 

~ [c 4 - rfic 2 - d ] - [c 4 - rfic 2 - d ] 



3 



(71) 



The constants 620, &21, &40, &6i and di, do are analogous to the constant a\ defined in 

equ. (59). We derive the relevant orthogonality relations in App. B. Finally, (Vf) zm is given 
by 



A/(°)(l + a/(°)) 



2T 



a\Tq l dia + 



26 



' 4 V 



15 



7d 2 i 
-7-9 -77^ 
6 3 



-a{Fq % diOi 



do, 



26 40 2 7d 2 .., 



+ 2c* 



ajdjP P 



(72) 



n n 

where we have used the Euler equations to eliminate time derivatives of the hydrodynamic 
variables. This is consistent at this order in the derivative expansion. The parameter 6 is 
defined in App. A (for a Maxwell gas 6 = n/T). We now use the Navier-Stokes equations 



= 2Ta + 2TV k d k a 



2b. 



a.Tq^a + =-^V 2 + - ^V^ 

15 6 3 



+ A 

= v i + v k V k Vi + Tdia 



'40 2 ^0 2 ^0, 



aiqi 







d t \ogT-V i q i + -(V-V+^- 



+ A 



9 



o\djP P^ k 

- LJ — + -v fe ^ 

n n 



(73) 
(74) 



—a{Fq % dia — 



2<? 

2&40 _2 ^0 „2 , ^0 

T5~ 



6 3 



to confirm that (Vf) zm + + = 0. Solving equ. (67) for /( 2 ) we get 



/( 2 ) = -A[(P/) orth + (P/) 



seal 



(75) 



(76) 



where (T>f) orth and {Vf) SC ai are defined in equ. (70) and (71). We note that (Vf) orth depends 
on Dlog(A). This term is necessary in order to respect conformal invariance. Dimensional 
analysis implies that A = T -1 /i(a:) for some undetermined function h. Using the Euler 
equation we can then show that Dlog(A) = |(cr). 
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D. Stress tensor and entropy current at second order 



Using the result for obtained in the previous section we can compute the dissipative 
correction to the stress tensor at second order in the gradient expansion. As discussed above 
only {Vf) or th contributes to the conserved currents. We will write jf\ h = —X(Vf) orth . We 

(2) 

can simplify yj th by combining some terms, and by using the Euler equation. We find 

2 



(2) A 2 /(°>(l + q/(°)) 

J orth ' 



c l c J c k c l a i] a kl - ^rcV 



+ (c 2 - a,) ((3c 2 - atf - 1) 



c 2 q 2 



+ [cWij - cVft + a^qi] -{a) + 2 (arf + 1 - f3c 2 ) c t c j c k a ij q k 



+ 
+ 



i j k j - fc 

CCC — C fyf 7 



VfcOij + [c 2 c fc - aic fe ] (q l V k Vi + £ ifc g' - g* - V'Vig*) 



cV 



c/^c 2 



(3n 2 - hP6)T 



+ 



-qidjd - (c 2 - ai) V^- 



(77) 



The stress tensor can be determined using equ. (56). The result is 



A 



V fc I — q k 
15 \ A 



(78) 



This result is in agreement with the general form of SU 1 ^ derived in Sect. HE. We note 
that some of the terms allowed by symmetry vanish in the kinetic theory calculation. This 
includes the Q* fe Q jfe structure and the VjPVjP and VjVjP terms. The kinetic theory also 
does not give a pure curvature term proportional to Rij. We also note that the coefficient 
of the g-> k a l k term can be written as trI], where tr is the relaxation time for the dissipative 
stresses. Our result shows that t r = A = rj/P. This result is in agreement with the 
calculation of spectral function of the stress tensor in [43], and with previous work on 
relaxation effects in dilute Bose and Fermi gases [44, 45]. 

Finally, we compute dissipative corrections to the entropy current and the dissipative 
function. We will use the expression for f s and R that are obtained by taking moments 
of the Boltzmann equation with log(//(l + a/)), see Sect. Ill A. As a consequence of the 
H-theorem these equations automatically satisfy the second law of thermodynamics. This is 
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no longer the case if the result is expanded in derivatives of the thermodynamic variables. 
The distribution function at second order in the gradient expansion determines R to third 
order in gradients, and this expression is not manifestly positive. It is possible, however, to 
add certain fourth order terms and obtain a manifestly positive result. This is the strategy 
we have adopted. We find 



3s 



( 27 Kq k a k * 2k 



Kq k tt ik Kg* (a) 



(79) 



and 



R 
T 



a- A 
+ kT (q- A 



1 2 

-aa -att + -(a) a + a + V k V k a 



27 qa qCl q(a) . k 



(80) 



In equ. (80) we have suppressed tensor indices. Contractions are defined as era = <7 ik cr k j, 
ail = a ik Vl k j etc., and (a) 2 = a^a^ . We note that the dissipative function is manifestly 
positive provided the first order relations r\ > and k > are satisfied, and that the 
entropy current is of the form discussed in Sect. II F. 



IV. OUTLOOK AND CONCLUSIONS 

In this work we have studied the constraints imposed by conformal invariance on the form 
of the hydrodynamic equations at second order in the gradient expansion. We find that the 
most general form of the stress tensor is given by 

5% = + r]T R (^gad) + V k V k a tJ + ^(a)a i3 ^j + Xia {i k a j)k + \ 2 v {i k tt j)k 

+ \ 3 n {t k n j)k + 7i v (i tv,)T + 72 v (i pv i> p + lz v { {rv j} p 

+ 7 4V (i V j> T + 7 5 V (i V,)P + K R R {lj) , (81) 

where we have suppressed terms involving the gauge fields. Conformal symmetry constrains 
the form of the comoving derivative of <Jij and eliminates possible trace terms like gij<J k \o kl . 
For phenomenological applications it is useful to rewrite the second order equations as a 
relaxation equation for the viscous stress 7Tjj = 5Uij. For this purpose we use the first order 
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relation 7r^ = —rjaij and rewrite equ. (81) as 

TCij = -rfffij - Tr ^ifc7T* + V fc Vfc7Tij + ^((j)^-^ + Ai7T (i fe 7r j)fe + . . . , (82) 

where . . . refers to the terms proportional to A2,3,7i and Kr. In deriving equ. (82) we have 
also used rj(n,T) = nf(a), where / is a function of a, and the Euler equations for n and 
a. Equ. (82) has the same structure as the second order equations considered by Israel and 
Stewart [29], but the coefficient | is specific to the non-relativistic conformal case. 

We have checked the conformal constraints by computing the dissipative contribution to 
the stress tensor, the entropy current, and the dissipative function in kinetic theory. We 
find that 511^ is indeed of the form given in equ. (81), but that some terms allowed by the 
symmetries do not appear in kinetic theory. The relaxation time for the dissipative stresses 
is given by t r = rj/P. The entropy current and the dissipative function also satisfy the 
conformal constraints. As in the case of the stress tensor, not all possible terms appear. 
The relaxation time for the entropy current is equal to the viscous relaxation time. 

There are a number of issues that we have not addressed in this paper. We have studied 
the constraints that arise from conformal symmetry and Galilean invariance, but we have 
not investigated the conditions that arise from the second law of thermodynamics. These 
conditions should have the form of inequalities for the transport coefficients, analogous to 
the first order relations rj > and k > 0. In the kinetic theory calculation the second law 
is satisfied automatically, but in the simple model considered here the transport coefficients 
are all governed by a single relaxation time A. We have also not attempted to derive Kubo 
relations for the coefficients tr, \, and 7,. The analogous Kubo formulas in the relativistic 
case were obtained in [32]. Finally, it would be interesting to compute the spectral functions 
for the stress tensor and the energy current at second order in the derivative expansion. 
These results would be useful in connection with attempts to compute transport coefficients 
using Quantum Monte Carlo methods [46]. 

Acknowledgments: This work was supported in parts by the US Department of Energy 
grant DE-FG02-03ER41260. We thank Dam Son for a useful discussion. 
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Appendix A: Thermodynamic relations 



In this appendix we collect a number of thermodynamic relations. In kinetic theory the 
number density and the pressure are determined by equ. (43) and (44). We have 



n(T, a) 



dXc/(0) (2T" a 



( c 2 



P(T,a) = \Jd Xc c 2 f^ (^-«) > 
where d Xc is defined below equ. (55) and we have set m — 1. We also define 



Q(T,a) 



dx c c 4 / (0) 



2T 



— a 



(Al) 
(A2) 

(A3) 



We can now compute various differentials. Differentials of n are given by 



dn(T, a) 



dT 
dn(T, a) 



2T 2 



dx c c 2 / (0) (l + a/^) = ^,(A4) 



{/dxJ (0) 

da {I dXcf(0) {^- a )} = l d Xcf {0 \l + af^) = 9T, (A5) 
where a = 0, ±1 corresponds to Maxwell, Bose-Einstein and Fermi-Dirac distribution 
functions, respectively. From the above relations we find 



3/7 

dn = 9T da + —d(\ogT) . 
Differentials of the pressure can be computed in the same way. We get 



dP(T, a) 



dT 
dP(T, a) 



-d T 
3 



da 



1 



2T 



( c 2 



5P 



nT, 



0„{ I A Xc c 2 r^--a 



and 



5P 3nT 15P 

dP = nTda + — d (logT) , d£ = —da + — — d (logT) . 



From £ + P = ytm + Ts we obtain 



ds(T, a) 



dT 
ds(T, a) 



5P 5 dP(T,a) 



da 



2nT 2 2nT dT 
5 <9P(T,a 



5P dn(T,a) 



-1 + 



2nT 9a 



2n 2 T dT 



5P dn(T,a) 



2n 2 T da 



3n 2 - 5P# 
2^2 



and 



3n 2 - 5P# , 
ds = da . 

2n 2 



For a Maxwell gas this relation becomes ds = —da. The differential of Q is given by 

70 

dQ = 15PTda + -^-d (logT) . 



(A6) 

(A7) 
(A8) 

(A9) 

(A10) 
(All) 

(A12) 
(A13) 
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Appendix B: Orthogonality relations 



In this appendix we collect the orthogonality relations that determine the coefficients 
b 2 o, b 2 i, 640, &41, &60> ^61 an d di, do introduced in Sec. Ill C. The coefficients b 20 , b 2 \ are defined 
by 

= / d Xc / (0) (l + «/ {0) ) ( i± ^^ 2 - ^ ~ «*>) , (Bl) 
= I d Xc / (0) (l + «/ (0) ) ( ^^ c 4 - & 2 ic 4 - 6 20 c 2 ) . (B2) 



These relations imply 



621 " 3n 2 -5Ptf' 620 " 2(3n 2 - 5P£) " (B3) 



The remaining coefficients are determined analogously. We have 




= / d Xc / (0) (1 + a/ (0) ) ( 1± ^c 4 - 6 41 c 2 - 6 40 ) , (B4) 
= J d Xc / (0) (l + «/ {0) ) ( 1+ 2 r /(0) c 6 " &4ic 4 - 6 40 c 2 ) , (B5) 



5(3n 2 - 7P#) 15Pn 



^ 41 " 2(3n 2 - 5P6) ' W4 ° " 3n2 - 5P6 ' (B6) 

= J d Xc + o /< >) ( 1+ 2 2 ° /(0) c 6 " &6ic 2 - 6 60 ) , (B7) 

= / d Xc /<°>(l + a/ {0) ) ( 1+ 2 ° /(0) c 8 - ^ic 4 - 6 60 c 2 ) , (B8) 

21(5nP-Qfl) 21(3nQ-25P 2 ) 

^ 661 " 2(3n 2 - 5P£) ' 660 " 2(3n 2 -5Ptf) ' (B9) 

J d X c/ (0) (l + o/ (0) ) (c 4 - rfic 2 - do) , (BIO) 

| d X c/ (0) (l + «/ (0) ) ( c6 - d ic 4 - d oC 2 ) , (BH) 

45nP - 7Q# 7nQ - 75P 2 

^ rfl " 3(3n 2 - 5Ptf) ' 3n 2 -5P£ " (B±2) 

where the thermodynamic variables 0, n, P, Q are defined in the previous section. 
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